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We address the in-plane pressure-dependent electrodynamics of graphite through synchrotron 
based infrared spectroscopy and ah initio Density Functional Theory calculations. The Drude term 
remarkably increases upon pressure application, as a consequence of an enhancement of both electron 
and hole charge densities. This is due to the growth of the band dispersion along the fcz direction 
between the K and H points of the Brillouin zone. On the other hand, the mid-infrared optical 
conductivity between 800 and 5000 cm~^ is almost flat, and very weakly pressure dependent, at 
least up to 7 GPa. This demonstrates a surprising robustness of the graphene-like universal quantum 
conductance of graphite, even when the interlayer distance is significantly reduced. 
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There is no doubt in considering the discovery of 
graphene [l| as one of the most important events in mate- 
rials science, during the last decade. The reason for this 
lies not only on the incredible number of potential ap- 
plications of this material, but also on the experimental 
observation of a completely new physics, which while al- 
ready being theoretically predicted 0, was not believed 
to be realized in any existing material. Graphene amazes 
with such a richness of phenomena occurring in one of the 
most simple existing crystal structures. 

Thanks to the discovery of graphene, also the research 
on its "real- world" analogue graphite, took momentum. 
Interestingly, many of the unusual phenomena discovered 
in graphene, such as the linear (Dirac-like) band disper- 
sion at the H point 0, the anomalous quantum Hall ef- 
fect Q , or the universal conductance properties , have 
been found to be present in graphite as well [6|- Un- 
derstanding the transition from graphene to graphite is 
therefore of the uttermost importance. It is within this 
framework that one should address the role played by the 
interlayer interaction in graphite, and ask what happens 
as long as the interlayer distance is externally modulated. 
The interlayer distance, for instance, can be expanded, 
by intercalating graphite with Br Q or K Q. Alter- 
natively one may try to study the behavior of graphite 
when the interlayer distance decreases and therefore the 
interaction among graphene layers enhances. 

This is indeed the strategy we chose, by addressing the 
pressure-dependent electrodynamics of graphite with in- 
frared reflectivity measurements and with new ah initio 
Density Functional Theory calculations. In particular, 
we aimed at noticing if pressure was able to demolish 
the universal quantum conductivity of graphite, previ- 
ously observed by Kuzmenko et al. Q . As in the case of 
graphene, the optical conductivity of graphite at ambi- 
ent pressure is ruled by transitions between hole and elec- 
tron bands, yielding a universal conductance Go = /Afi 



per layer, which is almost frequency independent be- 
tween 800 and 5000 cm^^. As we will see in the fol- 
lowing, pressure affects the out-of-plane band dispersion 
and the low-frequency conductivity of graphite by in- 
creasing the Drude spectral weight. Moreover, also at 
pressures as large as 7 GPa, a huge portion of the mid- 
infrared conductivity is practically independent of pres- 
sure. This shows that the universal quantum conduc- 
tance of graphite is robust against a strong increase in 
the interlayer interaction. 

A Highly Oriented Pyrolitic Graphite (HOPG) chunk 
(mosaic spread angle < 0.4°), purchased from SPI Sup- 
plies of about 100 micron size was loaded into a screw- 
driven Diamond Anvil Cell (DAG). The sample was po- 
sitioned in the DAG such as to obtain a flat and clean 
interface between the ah plane of graphite and diamond. 
Gsl salt was used as pressure transmitting medium, while 
a ruby chip loaded in the cell together with the sample is 
used as pressure calibrant through the fluorescence tech- 
nique. Measurements have been performed with the help 
of a Bruker 66v spectrometer and Hyperion 2000 infrared 
microscope over an extended frequency range, between 
350 and 13000 cm~^. This was made possible by the use 
of synchrotron radiation Q, as a bright source for far- 
and mid-infrared microspectroscopy. 

The ambient temperature reflectivity of graphite basal 
plane at the sample-diamond interface {Rsd{^^)) is re- 
ported in Fig. [1^. The break in the experimental 
data points between 1700 and 2500 cm"-'^, is due to 
diamond multi-phonon absorptions. Rsd{'-^) monoton- 
ically decreases with increasing frequency, except for 
a weak bump located at about 400 cm^^. With in- 
creasing pressure, the slope of i?sd(w) gets steeper in 
the spectral range below 3000 cm~^, thus suggesting 
an increased metallicity, while the bump gets more pro- 
nounced and shifts towards high frequencies. Reflectivity 
data were first fitted within a Lorcntz-Drude phenomeno- 
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FIG. 1: Color online, a) Graphite reflectivity in the far 
and mid-infrared range at sample diamond interface at 1 
and 7 GPa. Dashed lines show the extrapolated behav- 
ior, obtained from a Drude-Lorentz phenomenological fit, for 
Kramers-Kronig Transformations, b) Full range reflectivity 
at same pressures, c) Real part of the optical conductivity 
from Kramers-Kronig Transformations. The corresponding 
Drude-Lorentz fits (dashed lines), barely distinguishable from 
the Kramers-Kronig curves, are also plotted. Blue, green and 
pink curves correspond to the main components of the fit (see 
text). The black dotted line is a reference for the universal 
conductance Go expected for graphite (see Ref. 0]). d) Dif- 
ference of the optical conductivities between 1 and 7 GPa. 
e) Full range optical conductivity at 1 and 7 GPa, compared 
with ambient pressure data 0|. 



logical framework [T3, [HI (dashed lines in Fig. [T^), and 
then extrapolated in the missing ranges. This allowed 
to perform Kramers-Kronig Transformations, and calcu- 
late the real part of the optical conductivity ai{uj) ( Fig . 
[T];), by following a procedure described elsewhere 
The quality of the Drudc-Lorcntz fit can be verified in 
both panels a and c of Fig. [1] Fig. [1]: also allows to 
visualize the main fitting components: one Drude term 
(blue), one narrow component at about 400 cm^^ tak- 
ing into account the bump in reflectivity (green), one 
strongly overdamped oscillator mimicking the graphene- 
like universal conductance (pink). A fourth component 
(not shown here) has also been employed to take into 
account electronic transitions in the near- infrared 0, [l3| . 



As previously found by Kuzmenko et al. Q at am- 
bient conditions, the optical conductivity at 1 GPa of 
graphite is flat in a large portion of the infrared range, 
and its value per layer matches that of the frequency inde- 
pendent universal conductance Go = e^/4?i of graphene. 
Deviations from Go are due to the presence of a narrow 
Drude term, and to a peak found at 400 cm~^. On the 
low frequency side, the optical conductivity extrapolates 
at about 10000 (f2cm)~^, in fair agreement with dc data 
141 . With increasing pressure the Drude term increases 
as well: This is the direct consequence of the Rsdi^^) 
enhancement described before. Moreover, the interband 
conductivity peak at 400 cm~^ also increases and slightly 
blue-shifts. At larger frequencies, and up to the higher 
energy interband features (> 6000 cm~^) [15|, the be- 
havior of ai{Lu) is the same as for graphene, even at pres- 
sures as high as 7 GPa. From the point of view of the 
optical spectral weight (SW), the pressure-dependent in- 
crease of the low-energy (Drude-f peak) spectral weight is 
fully recovered at frequencies of the order of 8000 cm~^, 
as shown in Fig. [TJi, where the difference of the optical 
conductivities between 1 and 7 GPa is represented as a 
function of frequency. 

The most striking feature of our data is that while the 
overall graphite conductivity appears to be substantially 
unaffected by pressure, the Drude term enhancement is 
extremely significant. The in-plane Drude plasma fre- 
quency ujp, increases with pressure by a factor of 30% 
(from 4700 to 6100 cm"^), between 1 and 7 GPa (see 
red circles in Fig. [2^). Note that while we can assume 
a relatively large error bar of about 15% in the absolute 
estimate of Wp due to some arbitrariness in the Drude- 
Lorentz deconvolution, the relative error in the ujp pres- 
sure dependence is well below 5%. On the other hand we 
know from literature [l6| that between the same two pres- 
sure values, the intra-layer lattice parameter a decreases 
by 0.5% only (for comparison, the inter-layer lattice pa- 
rameter c decreases by 10% in the same pressure range). 
This rules out the possibility that a bandmass reduction 
due to the intra-layer compression may be responsible 
for the observed increase of the plasma frequency. The 
enhancement of Up has rather to be ascribed to a varia- 
tion in the number of the charge carriers upon pressure 
application. The microscopic mechanism underlying this 
effect will be investigated in the following, with the help 
of pressure-dependent ab initio calculations. 

The present ab initio calculations are performed in 
the local density approximation to the density-functional 
theory (LDA-DFT) as implemented in the all-clcctron 
LAPW code Wien2k [l^. A muffin tin radius of 1.3 Bohr 
is used, with a Rmt x Knax — 7 plane wave cutoff. Self- 
consistency is achieved by Brillouin Zone integration on 
a 60 X 60 X 20 Mokhorst-Pack mesh. We use experimen- 
tal lattice parameters according to Halfland et al.fl^. 
The use of augmented plane waves has the advantages 
that wavefunction norm conservation provides a direct 
access to the computation of the linear momentum oper- 
ator pij by the expectation value of the gradient operator 
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As outlined above, this effect can not be attributed to 
a bandwidth increase (i.e. a mass reduction) due to the 
shght in-plane compression. To validate this hypothe- 
sis we have performed ab initio calculation for a ficti- 
tious graphite system, subject to an in-plane compres- 
sion, but with the same intcrlayer spacing as ambient 
pressure graphite. The calculation (triangles in Fig. [5^) 
shows that under these conditions ujp docs not increase 
upon intra-layer compression, but even slightly decreases. 
This tells us that the intra-layer compression induced by 
pressure has only a negligible effect on the in-plane elec- 
tronic properties. 



FIG. 2: Color online, a) Pressure dependent plasma fre- 
quency from data (solid symbols) and calculations (open sym- 
bols) plotted as a function of pressure. The plasma frequency 
calculated in the case of pure intralayer compression (no vari- 
ation of interlayer parameters) is also shown (solid triangles), 
b) Electron (blue) and holes (red) FS sheets in ambient pres- 
sure graphite. 



{~ihV)ij. This is especially useful for the determination 
of the optical conductivity we compute as a sum of a in- 
traband (Drude) and interband term. The Drude's term 
is computed from the knowledge of the plasma frequen- 



cies ojp^a (a = x{y) or z) defined as: 



\Pv. 



(1) 



For the interband term we calculate the (complex) di- 
electric function deriving the optical conductivity from 
the imaginary component given as: 



«(w) = 



1 
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■^J.a\^5{£i{k)-£J{k)-nu:) , 

(2) 

where ej{k) < Ep < Si{k). The accuracy in the optical 
conductivity calculations comes from the correct descrip- 
tion of the Fermi Surface (FS) topology and convergence 
of the BZ summation in Eqs. ([T]) and ([2]). As for the 
FS topology we obtain for the ambient pressure graphite 
two sheets (electrons and holes) shown in Fig. |3)d. The 
shape for the electron (blue) and holes (red) is in agree- 
ment with the results of Refd andfiol The convergence 
of the BZ summation is very important in computing the 
plasma frequency. To get a confidence within a few per- 
cent we found necessary to go up to a mesh equivalent to 
a 600 X 600 x 100 Monkhorst-Pack. 

Knowledge of the complex dielectric function allows to 
recalculate the reficctivity at the sample diamond inter- 
face, at any given pressure. As reported in Fig. [31 it 
can be easily seen that the overall agreement between 
the calculations of Rsdi^^) and cri(w) with data is good. 
The plasma frequency values, both from data and cal- 
culation, are reported in Fig. [5^, as a function of pres- 
sure. It is immediately seen that the steep enhancement 
of ujp described previously, is confirmed by calculations. 
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FIG. 3: Color online. Comparison of the pressure dependent 
optical data with the results of the ab initio calculation in 
terms of reflectivity at sample diamond interface (a,b) and 
optical conductivity (c,d). 

In Fig. 2^ we show the band structure dispersion for 
free standing and compressed (7 GPa) graphite along the 
M-K'-H'-L directions, where K'-H' stands for a line par- 
allel to the K-H but slightly shifted towards the M point 
(see lower inset). This choice allows to see more clearly 
the bonding tt and antibonding tt* bands that are de- 
generate along the KH direction (see upper inset). The 
Figure shows the origin of the holes and electrons Fermi 
surfaces sheets of Fig. [^b from the tt and tt* bands re- 
spectively. Pressure broadens the bands along the k^ di- 
rection between points K' and H' in Fig. |4^, while leaving 
the in-plane dispersion (along K'M and H'L) almost un- 
changed. The increased kz dispersion moves the K and H 
semi-metallic points further apart in energy, passing from 
37 meV at ambient pressure to 112 meV at 7 GPa, with 
the chemical potential lying in between. This results at 
K and H, in a relative shift - of opposite sign - of the 
semi-metallic crossings with respect to the chemical po- 



4 



tential, thus leading to an increased number of free charge 
carriers of both signs (Fig. HJj). This quahtative assess- 
ment is confirmed by a more detailed calculation show- 
ing how the number of free electrons and holes increases 
with pressure by the same amount. Moreover compres- 
sion is also responsible for the shift of the peak around 
400 cm~^ in the optical conductivity. Indeed this peak 
results from the interband transitions shown in Fig. 2^ 
(vertical arrows) . Along the KH direction this absorption 
does vanish because of the degeneracy of tt and tt* bands. 
In the neighboroud of that direction (namely K'H') the 
bands are non degenerate and their contribution reaches 
its maximum value, thus giving rise to a peak in ai{u;), 
where tt and tt* are full and empty respectively. The 
direct gap along the K'H' direction strongly depends on 
interlayer distance and justifies the upshift of the peak 
in the optical conductivity as experimentally observed. 




FIG. 4: Color online. Graphite band structure dispersion at 
and 7 GPa from ab initio calculations, a) Band structure 
along the M-K'-H'-L directions (see text). The figure high- 
lights the mechanisms responsible for the pressure dependent 
increase of carrier concentration, as well as for the blue shift 
of the interband absorption peak experimentally observed at 
about 400 cm~^ (see text), b) Same along the high symme- 
try lines M-K-H-L, where the bonding tt and antibonding tt* 
bands get degenerate. 



As a final remark, we comment on the pressure- induced 
spectral weight redistribution from the far to the mid- 
infrared, as shown in Fig. [TJi. Interestingly, this behavior 
is q ualitatively similar to that observed by Horng et al. 
[20| in biased graphene. These authors find a transfer of 
spectral weight from low (< 1000 cm~^) to high (> 1500 
cm~^) frequencies. In graphene this behavior is explained 
by the electric-field tuning of the chemical potential close 
to the Dirac point: When the chemical potential moves 
away from the Dirac point, the Fermi Surface (and thus 
the Drude term) increases at the expense of the low en- 
ergy interband transitions @ . In a similar way, pressure 
in graphite, increases the number of charge carriers, while 
depleting the low energy interband transitions following 
the same mechanisms as described in Refs. |^ [13 and 
[21I . To this aim, the only band-structure element needed 
is the presence of a band crossing close to the chemical 
potential, independent from the details of the band dis- 
persion (either conic as in graphene or parabolic as in 
graphite K point). 

We have shown in this letter how pressure affects the 
in-plane low energy electrodynamics of graphite, by in- 
creasing the number of both electron and hole charge 
carriers. This effect is due to a pressure-induced shift of 
opposite sign, of the bands at the K and H point, with 
respect to the chemical potential. The other effect of 
pressure is the enhancement and blue shift of an inter- 
band transition located at about 400 cm~^. Above this 
frequency, the effects of pressure on the in-plane opti- 
cal conductivity are rather limited, at least up to 1 eV, 
where the inter-band transitions previously studied by 
Hanfiand et al. [l^ take place. Nonetheless, a very large 
portion of the mid-infrared spectral range between 800 
and 5000 cm~^ is very weakly pressure-dependent, thus 
demonstrating a surprising robustness of the universal 
conductance properties of graphite when the interlayer 
distance is significantly reduced. 
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